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ABSTRACT 

The probability density function (PDF) of the gas density in subsonic and supersonic, isothermal, 
driven turbulence is analysed with a systematic set of hydrodynamical grid simulations with resolutions 
up to 1024 3 cells. We performed a series of numerical experiments with root mean square (r.m.s.) Mach 
number M. ranging from the nearly incompressible, subsonic [Ai — 0.1) to the highly compressible, 
supersonic (AA = 15) regime. We study the influence of two extreme cases for the driving mechanism 
by applying a purely solcnoidal (divergence-free) and a purely compressive (curl-free) forcing field to 
drive the turbulence. We find that our measurements fit the linear relation between the r.m.s. Mach 
number and the standard deviation of the density distribution in a wide range of Mach numbers, 
where the proportionality constant depends on the type of the forcing. In addition, we propose a 
new linear relation between the standard deviation of the density distribution a p and the standard 
deviation of the velocity in compressible modes, i.e. the compressible component of the r.m.s. Mach 
number A^ C omp- In this relation the influence of the forcing is significantly reduced, suggesting a 
linear relation between <j p and A^ C omp, independent of the forcing, ranging from the subsonic to the 
supersonic regime. 

Subject headings: compressible turbulence 



1. INTRODUCTION 

Understanding the intricate interplay between in- 
terstellar turbulence and self-gravity is one of the 
key problems in star formation theory. The super- 
sonic turbulent velocity field is likely responsible for 
the complex and filamentary density structures ob- 
served in molecular clouds. It creates dense re- 
gions that can become gravitationally unstable and 
collapse i nto dense cores and eventually turn into 
new s t ars (lElmeereen fc Scaloll2004l : lMac Low fc Klessenl 
120041 : iMcKee fc Ostrikerl 120071 ). Statistical quanti- 
ties, describing this process, such as the initial mass 
funct i on, the core mass funct i on (IPadoan fe Nordlundl 



the mass-density grows proportional to the root mean 
square (r.m.s.) Mach number A4 of the turbulent flow 



20091) . and the 



l200l iHennebelle fc Chabrierl 

star formation rate flHennebelle fc Chabrierl 120111 : 
iPadoan fc Nordlundl 1201 1[ ) depend on the standard de- 
viation (std. dev.) of the density of the mole cular 
cloud. The pioneering works by IPadoan et ail (|1997l ) and 
iPassot fc Vazquez-Semade"m1 (|1998f ) have shown that the 
std. dev. a p of the probability density function(PDF) of 
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where (p) v is the volume-weig hted mean densi t y and 
6 is a proportionality constant. iFederrath et al.l (|2008l 
I2010T ) explained the dependence of a p on b by taking 
the modes of the forcing into account that drive the 
turbulent velocity field. This model predicts for purely 
solcnoidal forcing 6=1/3 and for purely compressive 
forcing 6=1, and explains the large deviations of 6 rang- 
ing from 6 = 0.26 t o 6 = 1.05 in previous works (e.; 



Padoan et al.l I1997L IPassot fc Vazquez-Semadenil I 199E 
Li et all 12003]: iKritsuk et all 120071: iBeetz etaU l20Qg 
Schmidt et"aLl 120091: iPrice et alJl2011b iKonstandin et all 
2012tlMolina" et al.ll2012f ). Wefollow up on this work and 



discuss the physical origin of this dependency and intro- 
duce a new relation, similar to equation ([1]), however, cor- 
relating the compressible component of the r.m.s Mach 
number A'Jcomp with a p . 

In section [2] wc explain our numerical setup. We analyse 
the influence of measuring mass-weighted and volume- 
weighted distributions in section |3~T1 the influence of the 
resolution on our measurements in section 13.21 and the 
PDFs of the mass density and the compressible part of 
the velocity field in section l3~3l In section l3~4"l we present 
the linear relations between the std. dev. of the mass den- 
sity and the r.m.s. Mach number. In section [3. 51 we dis- 
cuss the new relation between the std. dev. of the mass 
density and the std. dev. of the compressible part of the 
velocity field. A summary of our results and conclusions 
is given in section 01 

2. SIMULATIONS AND METHODS 

We use the piecewise parabolic method 
(|Colella fc Woodward! [I9841) implemen t ed in the 
grid code FLASH3 (jFrvxell et all 120001 ; iDubev et al] 
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l2008f ) to solve the hydrodynamical equations on a 
uniform three-dimensional grid. These equations are the 
continuity equation 



|^ + (vV)p = -pV-v, 



(2) 



the Euler equation with a stochastic forcing term F per 
unit mass 



— + v-V v = ^ + F, 

ot p 



and the equation of state 



p = up , 



(3) 



(4) 



where v is the velocity field, s = ln(p/ (p) v ) is the natu- 
ral logarithm of the mass density p, c s is the sound speed, 
p the pressure, T the adiabatic index. Since isothermal 
gas is assumed throughout this study, T = 1, the pres- 
sure, p = pc%, is proportional to the mass density with 
a fixed sound speed c s . These simulations are scale free, 
so we set (p)v = 1 5 c s = 1, and the box size of the sim- 
ulation L = 1. The numerical simulations are evolved 
for ten dynamical time scales T = Lj (2Aic s ), where 
•M = v r m. s Jc s is the r.m.s. Mach number of the simula- 
tions with the r.m.s. velocity v r . m .s.- All relevant quanti- 
ties are stored in intervals of 0.1T. The stochastic forcing 
field F has an autocorrelation time equal to the dynami- 
cal time scale on the injection scale, and varies smoothly 
in space and time. The forcing field is constructed in 
Fourier space such that the kinetic energy is injected on 
the largest scales, where 1 < kL/2-K < 3. To analyse 
the influence of different modes of the forcing field, we 
use projection tensors in Fourier space to get a purely 
divergence- free, V-F = 0, solcnoidal or a purely curl- free, 
VxF = 0, compressive vector field for the forcing. We 
adjust the amplitude of the forcing, such that we have 
r.m.s. Mach numbers Ai = 0.1, 0.5, 2, 5.5, 15 for both 
types of forcing in the stationary state of fully devel- 
oped turbulence. To investigate the effects of numerical 
viscosity, we study simulations at different resolutions, 
128 3 , 256 3 , 512 3 and 1024 3 . The parameters o f thes e 
simulations are described in iKonstandin et all (|2012t ). 
and a detailed descr iptio n of the forcing i s prese nted in 
iSchmidt et all p009l ) and iFederrath et all (f20l"ol ). 

3. RESULTS 

Figure [T] shows the time evolution of the r.m.s. Mach 
numbers Ai in all simulations. The fluid reaches the equi- 
librium state of fully developed turbulence after about 
two turbulent crossing times t fs 2T. We thus average 
all the following analyses for 2 < t/T. 

3.1. Volume-weighted and mass-weighted probability 
density functions 

It is well-known that the PDF of the logarithm 
of the mass density p(s) in a turbulent, isother- 
mal med ium is close to a Ga ussian distribution 

(see e.g. iVazquez-Semadeni 119941 : iPassot et al.l 1 1994 ; 
Padoan et all Il997t iKlessenl 120001: iKritsuk et al.l 12007 : 
8t IKonstandin et al.ll2012[) 
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Fig. 1. — R.m.s. Mach number of all simulations as a function 
of the dynamical time scale, calculated by averaging over all grid 
cells for both types of forcing. 



I Li et al.l (|2003f) showed with the assumption of a Gaus- 
sian, volume-weighted PDF of s that the mass-weighted 
PDF of s is also Gaussian with the same std. dev. and 
with a shifted mean value 
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Figure [2] shows the volume- and mass- weighted PDFs 
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Fig. 2. — Mass- weighted and volume- weighted PDFs of the loga- 
rithm of the mass density in the simulations with Ai = 5.5, 1024 3 
grid cells and both types of forcing. 



M 



(the volume- weighted PDF is shifted with (s) 
a s for a better comparison) for the simulation with 
Ai = 5.5 for both types of forcing. The PDFs are av- 
eraged over 81 time snapshots in the state of fully de- 
veloped stationary turbulence for t 2T and the error 
bars indicate the std. dev. of the temporal fluctuations. 
The variance of the volume- weighted PDFs is larger than 
the variance of the mass-weighted distributions. This ef- 
fect is stronger for the compressive forcing than for the 
solcnoidal forcing. The volume-weighted PDFs show a 
larger variation with time in the low-density wing of the 
distribution than the mass-weighted distributions. This 
low-density wing also shows higher probabilities than one 
would expect from the underlying Gaussian distribution 
extrapolated from the high density wing. This effect is 
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stronger for the compressive than for the solenoidal forc- 
ing. We assume that this behaviour is caused by our forc- 
ing scheme. As the time correlation of the forcing field 
is equal to the dynamic time scale on the largest scales, 
the forcing has enough time to produce very low densities 
in large regions of diverging flows. This process causes 
the volume-weighted PDF of s to have a tail at low den- 
sities with higher probabilities than the distribution for 
the case of turbulence, which is not driven on the largest 
scales. With increasing r.m.s. Mach number and increas- 
ing amplitude of the forcing field this effect becomes more 
important and the low-density tail influences the calcu- 
lated std. dev. of these distributions. This effect is less 
pronounced measuring the mass- weighted distributions, 
as the very low density grid cells carry only little mass. 
We note that there are other potential processes, which 
could lead to non-Gaussian w ings in the P DF, such as 
turbulent intermittence (e. g. [Klessen 2000h. 
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mass density, it can be shown that th e std. dev. of the 
Gaussian-distributed quantity s is (see iPrice "elall[20rl 
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Fig. 3. — Mass-weighted PDFs of s of the simulations for A4 = 
5.5, different resolutions and both types of forcing. The black solid 
lines are Gaussian functions with mean value and std. dev. calcu- 
lated with the highest resolution. 



3.2. Resolution effects on the probability density 
functions 

Figure [3] shows the mass- weighted PDF of the quantity 
s with an r.m.s. Mach number Ai = 5.5 and different 
resolutions. The PDF of s shows deviations from the 
Gaussian shape and a dependency on the resolution only 
in the high-density tails of the distribution. We interpret 
the deviations of our measured PDFs from the Gaussian 
distribution in the supersonic regime for both types of 
forcing as a sign of numerical dissipation and finite sam- 
pling. In the highly supersonic regime the medium is 
dominated by shock fronts and high-density gradients, 
which require high resolution to converge. As we do not 
fully resolve them in the Ai = 5.5 case, an additional 
dissipation occurs. This effect is stronger in the simu- 
lations with compressive forcing and becomes stronger 
with increasing r.m.s. Mach number for both types of 
forcing (not shown here). However, increasing the resolu- 
tion has no influence on the deviations from the Gaussian 
distribution in the low-density tail of the mass-weighted 
PDFs. 

With the assumption of a log-normally distributed 



(7) 



Figure [4] shows o~ p as a function of a s for our volume- 
weighted (left panel) and mass-weighted (right panel) 
distributions. The volume- and the mass-weighted mea- 
surements of the std. dev. of s show increasing deviations 
from equation ([7]) with increasing r.m.s. Mach numbers 
for both types of forcing. However, the deviations are 
smaller in the mass-weighted case than in the volume- 
weighted one. The assumption of Gaussianity, which is 
implied in equation ([7]), is better fulfilled for the mass- 
weighted case. Figure U] also shows that our measure- 
ments with Ai = 15 are not converged with resolution for 
both types of forcing. Our m easurements are in agree- 
ment with IPrice etaLl (|2011l) , who showed that direct 
measurements of o~ p show a stronger dependency on res- 
olution than measurements of a s . 

All volume-weighted measurements show a clear trend 
towards the relation ([JJ with increasing resolution. How- 
ever, the data points do not fit relation (JTJ) for Ai = 15 
with solenoidal forcing and in all the supersonic cases 
with compressive forcing, although the data points with 
Ai = 2 and Ai = 5.5 with compressive forcing are nearly 
converged with resolution. This can be explained with 
the non-Gaussian, low-density wing of the distributions, 
which docs not depend on the resolution. Considering 
that the std. dev. ct S; a/ of the mass- weighted PDF is more 
compatible with the scaling for a log-normal PDF, equa- 
tion (UJ, and that the resolution dependence of o~ s ,ai is 
weaker that for (T s y, we prefer to use <7 Si m as estimate 
for the turbulent density fluctuations in the following. 

3.3. The probability density function of the density and 
of the compressible modes in the velocity field 

Figure [5] shows the mass-weighted PDFs of the quan- 
tity s (left panels) and the volume-weighted PDFs of 
the compressible modes of the velocity field normalised 
to the sound speed M comp = v conip /c s (right panels) 
for different r.m.s. Mach numbers and both types of 
forcing. The PDFs of the logarithm of the density 
largely follow Gaussian distributions for all supersonic 
r.m.s. Mach numbers. We added Gaussian functions 
(black solid lines), with the first- and second-order 
moments calculated from our distributions in figure 
The high-density tails of the distributions show 
deviations from the Gaussian shape, which increase with 
increasing r.m.s. Mach number. Also the deviations 
from the Gaussian distribution in the low-density tail, 
as discussed in section 13.11 get more pronounced with 
increasing r.m.s. Mach number. Thereby, we have 
large deviations of our measurement from the Gaussian 
distributions in the Ai = 15 case and the calculated 
std. dev. does not correspond to the std. dev. of the 
underlying Gaussian distribution. 

The density distributions of the simulations driven by 
solenoidal forcing in the subsonic regime show significant 
deviations from the log-normal shape, which become 
stronger as Ai decreases. These distributions become 
more asymmetric and more peaked. The different be- 
haviour of the PDFs in the subsonic regime especially for 
the solenoidal forcing is caused by the different physical 
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Fig. 4. — Std. dev. of the mass-density a p as a function of the std. dev. of the logarithm of the mass-density <r s , measured volume-weighted 
(left panel) and mass- weighted (right panel). The deviations of the measurements from the black solid lines, equation J7}, quantify the 
deviations from a log-normally distributed mass density. 

processes acting here. In the subsonic regime sound 
waves transfer information faster than the averaged flow 
of the medium, such that the thermal pressure increases 
before two converging flows can collide. This process 
prevents colliding flows from producing high-density 
regions and causes the sharp edge at the high-density 
wing of the distributions. The thermal pressure also 
decelerates the velocities in compressible modes, such 



that the PDF of M comp also shows a narrow, peaky and 
intermittent behaviour for the solenoidal forcing. This 
process is just visible for solenoidal forcing, because 
in the compressive forcing case the velocities in com- 
pressible modes are re-injected by the forcing to hold 
the r.m.s. Mach number constant. This is the reason 
why the thermal pressure does not have such a strong 
influence there. 

The right panels of figure [5] show the PDFs of M comp , 
where M comp is calculated in Fourier space and averaged 
over the three directions of the coordinate system x, y, z 
in real space. The distributions of Af comp are symmetric 
with zero mean and have an increasing std. dev. with 
increasing r.m.s. Mach number. The distributions 
obtained with compressive forcing are always broader 
than with solenoidal forcing at the same r.m.s. Mach 
number. The PDFs of M comp are Gaussian (black solid 
lines) with deviations in both wings. These are the 
signpost of turbulent intcrmittency. The deviations do 
not show a clear trend with the r.m.s. Mach number. 
The PDF of M comp obtained with solenoidal forcing in 
the subsonic regime with Ai — 0.1 shows the strongest 
deviations from the Gaussian shape with a narrow, 
peaky, intermittent distribution. These deviations are 
caused by the thermal pressure, as discussed above. 



3.4. 



In 



Relation between the r.m.s. Mach number and the 
standard deviation of the density 

iPadoan et al.l (|1997t ) and 



number dependency of density fluctuations is used to 
derive analytic expressions for the core mass function 
(CM F) and the stellar initial mass function (IMF) 
(e.g.. IPadoan fc Nordlundl[200llHennebelle k Chabrierl 
I2008L 120091 ). On galactic scales it i s used to rep roduce 
the Kenni c utt- Sc hmidt relation (jTassis! l2007t ). and 
lElmegreenl (|2008t ) suggests that the star formation 
efficiency is a function of the density PDF. Figure [5] 
(upper left panel) shows the measured std. dev. of the 
mass density as a function of the r.m.s. Mach number for 
different resolutions and both types of forcing. The error 
bars in each panel of figure [6] indicate the std. dev. of the 
temporal fluctuations of the measured quantities. They 
do not include any potential systematic errors stemming 
from, e.g., the numerical scheme or implementation of 
the forcing algorithm. Thus, we interpret the error bars 
as a lower limit of the real uncertainty. The dotted 
and dashed-dot t ed lin es correspond to the model of 
iFederrath et al.l (|2010D . which describes the proportion- 
ality parameter & as a function of the turbulent forcing. 
This model predicts for solenoidal forcing b = 1/3 
and for compressive forci ng b = 1. Our measu rements 
agree with the model of IFederrath et all (j2010l ) in the 
supersonic case for both types of forcing. We see small 
deviations from the model in the simulations with 
A4 = 15, which is caused by our limited resolution (see 
figure |4]). The std. dev.s of the density distribution of 
the simulation with solenoidal forcing are smaller than 
the prediction of the model in the subsonic case. In 
the subsonic regime, the deviations arc caused by the 
thermal pressure, which damps density variations and 
compressible modes of the velocity field and reduces 
the measured std. dev. below the model prediction as 
discussed in section 13.31 The upper right panel of figure 
[6] shows the mass- weighted, logarithmic std. dev. <t Sj m 
as a function of the r.m.s. Mach number. The dotted 
and dashed-dotted lines correspond to the standard 
model for the logarithmic density variance, 



iPassot k Vazquez-Semadenl (|1998| ) the authors found 
that the std. dev. of the PDF of the mass density a p is 
proportional to the r.m.s. Mach number in a turbulent 
flow. The std. dev. of the mass density is an important 
quantity especially in astrophysics, where the Mach 



a\ = ln(l + b 2 M 2 ) 



(8) 



with 6=1/3 for solenoidal and b = 1 for compressive 
forcing. Equation ([5J follows from equatio ns (Q and (Ht 
and was recently derived analytically by I Molina et al.l 
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Fig. 5. — The mass- weighted PDFs of the logarithm of the mass density (left panels) and the compressible part of the local Mach number 
(right panels) for different r.m.s. Mach numbers, resolutions and both types of forcing. In the inset, a magnification of the PDFs obtained 
with solenoidal forcing for M = 0.1 are shown. The error bars in each panel indicate the std. dev. of the temporal fluctuations. 
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Fig. 6. — The std. dev. of the distribution of the mass density (left) and the std. dev. of the distribution of the logarithm of the mass 
density (right) as a function of the r.m.s. Mach number (upper panels) and a s a function of the com pressible part of the r.m.s. Mach 
number (lower panels). In the upper panels the lines correspond to the model of Fcdcrrath ct al. (201(3) with 6=1/3 for solcnoidal forcing 
and 6 = 1 for compressive forcing. In the lower panels the solid lines correspond to a two-parameter fit and the dotted line corresponds 
to a linear relation between the std. dev. of the mass density and the std. dev. of the compressible part of the r.m.s. Mach number with a 
proportionality constant \/3. 

3.5. Physical origin of density fluctuations in turbulent 

flows 

Studying the continuity equation ([2]), one can argue 
that variations of the density can only be caused by the 
divergence of the velocity field. Given that a vector field 
can be decomposed in a gradient field and a rotation field 
and that the divergence of a rotation field vanishes, we 
conclude that the density variations can only be caused 
by the compressible modes of the velocity. For this reason 
we replaced the r.m.s. Mach number, which is in fact the 
std. dev. of the velocity distribution, with the compress- 
ible part of the r.m.s. Mach number, .Mcomp, hi the lower 
panels of figure |6l The data points show a clear correla- 
tion. The different behaviour of the simulations driven 
with solcnoidal and compressive forcing arc significantly 
reduced. We added in figure [6] a function (dotted line) 



((2011) using the shock-jump conditions and averaging 
over an ensemble of shock waves. The deviations of our 
numerical data from this standard model are only signif- 
icant for solenoidal forcing in the subsonic regime, while 
our data are in excellent agreement with equation ([5]) 
for both solcnoidal and compressive forcing in the su- 
personic regime, given our resolution dependence of the 
M. = 15 data po i nts. Our results arc in agreement 
with IKowal et al.l (|2007t) . who found deviations from 
the linear relation with a p in the subsonic regime with 
soleno idal forcing, and with lPassot fc Vazquez-Semadenil 
(|1998f ). who analysed one-dimensional simulations with 
only compressive forcing and 0.5 < M < 3 and found 
the linear relation between M. and a p with 6=1 



I Price et al.l (|2011|) analysed three dimensional simula- 
tions with purely solenoidal forcing and r.m.s. Mach 
numbers in between 2 < M. < 20 and found 6=1/3 
in excellent agreement with our result. As they did 
not analyse the subsonic regime with solcnoidal forcing 
they did not observe the large deviations in the sub- 
sonic regime. Our analysis complements these studies 
with measurements in both the subsonic and supersonic 
regime and for purely compressive forcing. 



for the relation a 



V3M C 



which is the simplest 



model for this relation assuming isotropy. The factor of 
y/3 is due to the fact that we use the distribution of the 
compressible modes of the velocity field averaged over 
the three directions of the coordinate system 
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(9) 
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= \/37Wcomp • 

A similar mo del has also been suggested by 

iFederrath et al.l (|2010f ) . where the parameter b in 

equations ([T]) and ((HJ was approximated by \[D (vp) 
with the spatial dimension D of the turbulence (D = 3 
in our case) and the so-called compressive ratio 
which is the ratio of compressible to total velocity 
fluctuations. This and our simple model fits the data, 
but shows deviations for the simulations with solenoidal 
forcing and the lowest and highest Mach numbers. The 
deviations for the M. = 15 simulation are again caused 
by the resolution dependency of a p . Additionally, we 
perform a fit of our data (black solid line) with two free 
parameters, <r p — a\/3A4^ omp for the density relation. 
We obtain a normalisation a = 1.0 ± 0.1 and a slope 
/3 = 0.85 ± 0.04. For the s-relation we transform the 
fitted function with equation ([7]). The measurements of 
the std. dev. of the density have larger deviations from 
the model as the measurements of the std. dev. of s. 
However, the model fits the measurements in both cases 
and provides a good description for the data points 
in the subsonic regime with solenoidal forcing, which 
are strongly influenced by sound waves. We conclude 
that the thermal pressure damps the velocities in 
compressible modes in a way that the relation between 
the velocities in compressible modes and the density 
variations in a turbulent medium is in a statistical equi- 
librium state, even if the medium is strongly influenced 
by sound waves. The deviation of the scaling exponent 
from the simple model can be interpreted as additional 
dissipative effects, which are proportional to A^comp- 
An example for these physical processes are shocks, 
which cause deviations from the log-normal distribution 
of the density PDF. However, systematic errors with a 
dependency on the r.m.s. Mach number could also cause 
deviations from the linear scaling and would be another 
possible interpretation for our fitted scaling exponent. 
The shown relation between the std. dev. of the density 
and the compressible part of the r.m.s. Mach number 
in principle enables us to measure the kinetic energy in 
compressible modes in giant molecular clouds, without 
knowing the absolute r.m.s. Mach number, the driving 
mechanism or the sound speed. The relations shown 
in the bottom panels of figure |6] are valid in both, the 
subsonic and supersonic regime. 

4. SUMMARY AND CONCLUSIONS 

We have investigated the influence of solenoidal 
(divergence-free) and compressive (curl-free) forcing on 
the PDF of the mass density in subsonic and super- 
sonic turbulence with a set of three-dimensional nu- 
merical simulations. We analysed the relation between 
the std. dev. of the mass density distribution and the 
r.m.s. Mach number. We found a new relation between 
the std. dev. of the mass density and the std. dev. of the 
compressible part of the velocity field. Our main results 
are as follows: 

• Compressive forcing yields mass density PDFs with 
std. dev.s proportional to the r.m.s. Mach number 
with 6=1. For solenoidal forcing, we measure 
b = 1/3 in the supersonic regime. Our findings are 



in agreement with previous studies, which however 
only explored different subsets of the full parameter 
space investigated here. We also found deviations 
of our measurements from the linear relation with 
solenoidal forcing in the subsonic regime. These de- 
viations from the linear relation can be explained 
with sound waves, which damp the faint compress- 
ible velocities and prevent the medium from pro- 
ducing over-densities. 

• We found a unique relation between the std. dev. of 
the mass density and the compressible modes of the 
velocity field with a fit to our data. Our new rela- 
tion is independent of the driving mechanism and 
still holds in the subsonic regime, where the flow 
is mainly influenced by sound waves. It does not 
show a strong influence on the resolution and other 
effects, which may cause a non-Gaussian distribu- 
tion of the density. 

• Our relation enables us for the first time to mea- 
sure the kinetic energy in compressible modes in 
units of the sound speed, without knowing the 
r.m.s. Mach number, the driving mechanism or 
the sound speed of the medium. This measure- 
ment can be used to distinguish between sub- 
sonic and supersonic compressive turbulent mo- 
tions. It will in principle allow us to measure 
the composition of the kinetic energy in the in- 
terstellar medium by combining independent mea- 
surements of the tot al r.m.s. Mach number (e.g. 
iBurkhart et al.l I2009T) and the std. dev. of the den- 
sity distribution dBrunt et al.l [20101 iBruntl l2010t 
iSchneider et al.ll2012D . 
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